Assisted hopping and interaction effects in impurity models. 
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We study, using Numerical Renormalization Group methods, the generalization of the Anderson 
impurity model where the hopping depends on the filling of the impurity. We show that the model, 
for sufficiently large values of the assisted hopping term, shows a regime where local pairing cor- 
relations are enhanced. These correlations involve pairs fluctuating between on site and nearest 
neighbor positions. 
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I. INTRODUCTION. 

The formation of local moments is the simplest 
and most extensively studied manifestation of strong 
electron-electron repulsion in impurities in metal a 1 : 2 :? . 
This phenomenon can also be observed in mesoscopic 
systems coupled to metallic leads 4 ^. In the opposite 
limit, when the "impurity" becomes a metallic grain 
whose electronic structure needs to be described by many 
narrowly spaced levels, the electron-electron repulsion 
leads to Coulomb blockade 6 . The intermediate regime 
when the energy scales of interest, the temperature, the 
strength of the Coulomb repulsion and the level spacing 
are not too different in magnitude is not so well under- 
stood. A similar situation arises in the study of strongly 
correlated systems. The analogous of the Anderson im- 
purity model is the Hubbard modetf^ for extended sys- 
tems, which leads to a Mott transition and the forma- 
tion of local moments at half filling. The opposite limit, 
when many levels within the unit cell have to be included, 
leads, generically, to ordinary Fermi liquid behavior. 

The simplest modifications of the extended Hubbard 
model taking into account the multiplicity of levels in the 
ions in the unit cell lead to assisted hopping term a 9 : 10 , 
which modify substantially the phase diagram, and tend 
to favor superconductivity. Analogously, if the internal 
degrees of freedom of the metallic grain are included, 
beyond the constant interaction term which leads to 
Coulomb blockade, one finds non equilibrium effectsii* 1 ^ 
which suppress the Coulomb blockade. These effects are 
similar to the formation of an excitonic resonance found 
in the excitation of core electrons^. One can also ex- 
tend the Anderson impurity model to include the effects 
associated to the existence of many orbitals. When the 
influence of these orbitals is described in terms of effec- 
tive interactions within the restricted Hilbert space of the 
usual Anderson model, one indeed finds assisted hopping 
terms 14 . The resulting model has been analyzed using a 
mean field, BCS like decoupling of the interaction term!*, 
and by the flow equation method 1 ^ (for a description of 



the method, se o 16 : 17 ), which is well suited to the anal- 
ysis of impurity modelsi^i^^i^. These works show 
a tendency towards local pairing away from half filling, 
in agreement with the studies of bulk system a 9 i 10 . We 
present here more accurate calculations, using the Nu- 
merical Renormalization Group Method^,, which char- 
acterize the low energy properties of the model in a nu- 
merically exact way. 



We examine the effect of the assisted hopping term on 
the local density of states and find that with increasing 
assisted hopping amplitude the peak characteristic to the 
mixed valence regime of the Anderson model gets broad- 
ened and shifted to negative frequencies. This result is 
consistent with the fact that the assisted term enhances 
the occupation of the local level, also found in the cal- 
culations. To demonstrate the effect of the new term on 
the formation of the local moment, we compute the local 
spin spectral function and find that the spin susceptibil- 
ity gets suppressed as the assisted hopping term is turned 
on. This suppression might be interpreted as the trace of 
the pairing. We also check the pairing correlations in a 
more direct way and calculate the related spectral func- 
tions. Our main result is that -although the effect of the 
assisted term on the pairing correlations on the local level 
could be explained qualitatively by just considering the 
renormalization of the level- the so-called off diagonal 
correlations (pairing between the (i-electron and the con- 
duction electrons located at the impurity position) are 
also enhanced. 



The paper is organized as follows: The model, and the 
method of calculation are described in the next section. 
Then, the main results are presented. Section IV gives 
the main conclusions of our work. 
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II. THE MODEL AND METHOD OF 
CALCULATION. 

A. The model 

We describe the effects of the internai degrees of free- 
dom of the impurity in terms of effective interactions 
defined within the restricted Hilbert space of the stan- 
dard Anderson impurity 1 ^. Using perturbation theory, 
this approach is justified when the typical level spacing 
within the impurity, A, is much smaller than value of the 
Coulomb interaction, U. Then, it is easy to show that 
the Anderson model is recovered in the limit A/U — > oo. 
We will study the Hamiltonian: 
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where Ha = Hk + Hi mp + Hh y b is the Anderson Hamil- 
tonian, and the assisted hopping terms are included 
in Hassisted- We have also defined n d<s — d\d s and 
n d = n d . + + rid.-. The parameter dV determines the 
strength of the assisted hopping effects. For the sake of 
simplicity we consider a flat band with constant density 
of states for the conduction electrons, with half band- 
width D. In the following we use D as the energy unit. 



B. The Numerical Renormalization Group 
Method. 



To compute different quantities numerically we use 
Wilson's Numerical Renormalization Group method 23,24 . 
In this method -after the logarithmic discretization of 
the conduction band - one maps the original impurity 
Hamiltonian onto a semi-infinite chain with the impu- 
rity at the end. One can show that as a consequence 
of the logarithmic discretization the hopping along the 
chain decreases exponentially, t n ~ A - ™/ 2 where A is the 
discretization parameter, and n is the index of the site 
in the chain. The separation of the energy scales due to 
the decreasing hopping provides the possibility to diag- 
onalize the chain Hamiltonian iteratively and keep the 
states with the lowest lying energy eigenvalues as most 
relevant ones. Since we know the energy eigenstates and 
eigenvalues, we can calculate thermodynamical and dy- 
namical quantities directly (e.g. spectral functions using 
their Lehman representation). 
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FIG. 1: Value of n<j as function of dV for V = 0.2, ed = and 
U = 0.05. 



III. RESULTS 

The model has four parameters, e d , U, V and dV. Pre- 
vious studiesi 4 * 1 ^ suggest that the most significant de- 
viations from the standard Anderson impurity model 
take place when the filling of the d orbital is such that 
1.2 < (n d ) < 1.8. We focus our attention to that regime 
as well. The filling is determined by e d , U and also by 
dV. The mean field decoupling of the assisted hopping 
term shows that the latter term contributes to the renor- 
malizationof the position of the level. The value of n d 
as function of dV is given in Fig.pQ. Since the assisted 
hopping term renormalizes the impurity level downwards, 
that term favors the occupancy of the impurity level. 

The fact that the assisted hopping term, dV, and the 
electron-electron repulsion U have opposite effects can be 
appreciated in Fig. [2] where the influence of the value of 
U on n d , for a fixed value of dV , is shown. 

Fig- El shows the imaginary part of theone particle 
Green's function associated to the localized level, defined 
as: 

-G d (u>) = V|(0|4H 2 5{uj-e m ) 

7T * — ' I 

m,cr 

+ ]r|<0K|m)| 2 <5(a; + e m ) , (2) 

where |0) (|m)) is the ground (mth excited) state of the 
full Hamiltonian with energy eo = (e m ). 

In presence of no assisted hopping, the local density 
of states exhibits a broad peak at the Fermi level which 
is characteristic to the mixed valence regime of the An- 
derson impurity model. As dV is increased from zero, 
the results show a broadening of the resonance, also con- 
sistent with the fact that, at the mean field level, the 
assisted hopping term modifies the effective one particle 
hopping. As one can see in Fig [[3], the resonance gets not 



3 




FIG. 2: Value of rid as function of U for V = 0.2, e<j = and 
dV = 0.15. 
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FIG. 3: Imaginary part of the impurity one electron Green's 
function, 1/tt G d (w), for V = 0.2, e d = and U = 0.05. 



only broadened but also shifted down to negative ener- 
gies. This fact could already be foreseen from the result 
for the occupation, since the occupation (at T — 0) is just 
the integral of the density of states for negative frequen- 
cies. This enhancement of the occupation must leave its 
trace on the local moment formation as well. 

The influence of the assisted hopping term on the for- 
mation of a local moment can be inferred from Fig. [3], 



which displays the spin-spin correlation function: 

Im q z s (u) = - |(0|n dT - n dl \n)\ 2 5{uj ~ e„) (3) 



This function, for dV = shows a pronounced peak at 
low energies, which tend to the Kondo temperature be- 
low which the local moment is quenched as U increases. 
At higher energies, Im Qg(u>) ~ w -1 , which corresponds 
to a Curie susceptibility as function of temperature. This 
1.5 
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FIG. 4: Spin-spin correlation function for V - 
U = 0.05 and different values of dV. 



0.2, e d = and 



behavior crosses over to a Im Q z g[oS) ~ to regime which is 
characteristic to the compensation of the moment. (Al- 
ternatively, one could say that for very high frequencies 
Im Qg(uj) ~ co^ 1 as a consequence of the constant time 
correlator of the spin for very short times, while for small 
frequencies Im g z s (to) shows a linear dependence on u> cor- 
responding to the ~ t~ 2 asymptotics of the correlation. 
As dV is turned on, this peak becomes broader and it is 
shifted towards higher energies indicating the suppression 
of the local moment. This suppression is also consistent 
with the previous results and might be considered as a 
fingerprint of the local pairing. 

To get a more direct insight to the relation of the as- 
sisted hopping and the pairing, we now analyze the possi- 
ble existence of pairing correlations in the model. Within 
the reduced Hilbert space of the Anderson's impurity 
model, we can define two types of such correlations: 



F(w) = X)k°K d Il n >| <5(w-e n ) + |(0|d T d i |n)| 2 <5(w + e„) 

n 

F'(cj) = J2{0\d\d\\n)(n\ (d T q - d;c T )|0}5(a; - e n ) + <0|d T <f;|n)<n| (d\c\ - d[c\) \0)6{u + e„) (4) 
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FIG. 5: Local pair correlation, F(cu) in Eq.(0J as function of 
dV for e = -0.025, U = 0.05 and V = 0.2. 
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FIG. 6: Local pair correlation, F(u>) in Eq.(0J as function of 
e for dV = 0, U = 0.05 and V = 0.2. 
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FIG. 7: On site pair susceptibility, T, as defined in Eq.@ 
as function of e for dV = (top), and as function of dV for 
e = -0.025 (bottom). 
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where c a = ^ fc c ka , representing the metal orbital closest 
to the impurity. The function F(u>) in Eq.QJ gives the 
magnitude of the on-site pairing, as present, for instance, 
in the negative U Anderson's impurity model. The func- 
tion F'(uS) describes off diagonal pairing. Virtual Cooper 
pairs resonate between the on-site position and that in 
which one component of the pair is at the impurity and 
the other is in the metal. 

We first analyze the function F(u>), which is a mea- 
sure of the on-site pairing. The function F(u>) is shown 
in Fig. [5] as function of the value of dV. For compari- 
son, we plot the same function for dV = for different 
values of e in Fig.[S]. The two functions show similar 



FIG. 8: Off diagonal pair susceptibility, T\ as defined in 
Eq.@ as function of e for dV — (top), and as function of 
dV for e = -0.025 (bottom). 



behavior, suggesting that the main effect of the assisted 
hopping term on the on-site pairing correlations is the 
renormalization of the level, since its effect is difficult to 
be distinguished from the effect of a change in the impu- 
rity electron level. 

We estimate the tendency towards pairing by using 
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Eqs.0} to define generalized susceptibilities: 



T = I dto 
lo 



T 1 = 



du) 



F{uj) + F(-lo) 
F'{uj) + F'(-uj) 



(5) 



These functions measure the tendency towards the dif- 
ferent types of pairing. Fig.[7j gives T as function of e 
for dV = 0, and as function of dV for e = —0.025. As 
it can be appreciated, the inclusion of the assisted hop- 
ping term dV does not change significantly the value of 
this susceptibility. Moreover, its magnitude is not too 
different from the value obtained for dV — 0. On the 
other hand, Fig.jH] gives T 1 as function of e for dV = 0, 
and as function of dV for e = —0.025. In this case, the 
effects of dV are quite significant. While of dV — this 
susceptibility remains close to zero for the whole range 
of values os e studied, it raises swiftly as soon as dV is 
finite. Hence, the assisted hopping term induces an off 
diagonal pairing term of the type described by T' given 
in Eq.lQJ. 

This is our main result. In contrast to the local pairing 
on the d-level -which can be understood as the conse- 
quence of the primary effect, the renormalization of the 
level- the off diagonal pairing correlations are increased 
by the enhanced hopping rate which allows stronger fluc- 
tuations if the spin of the d electron and the local conduc- 
tion electron have opposite orientations. There is another 
interesting feature of this result: Though the formation 
of the local moment is strongly suppressed and the im- 
purity has no well defined spin state, the off diagonal 
correlations are still enhanced. This would indicate that 



this type of pairing manifests itself in a dynamical way: 
The spin state of both the local level and the conduction 
electron state fluctuate, but in a correlated way. 



IV. CONCLUSIONS. 

We have studied, using the Numerical Renormalization 
Group method, the Anderson impurity model with a hop- 
ping which depends on the charge state of the impurity. 
This is the simplest model which includes information 
about the internal structure of the impurity, beyond a 
single, rigid, electronic state. 

Our results indicate that, when the assisted hopping 
is sufficiently strong, the model shows a crossover to a 
phase with off diagonal pairing correlations. It would 
be interesting to know if this regime could be realized 
in mesoscopic devices, where the internal structure of a 
quantum dot can be important. 
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